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ABSTRACT 

Sgr A* is a source of strongly variable emission in several energy bands. It 
is generally agreed that this emission comes from the material surrounding 
the black hole which is either falling in or flowing out. The activity must be 
CLi driven by accretion but the character of accretion flow in this object is an 

JL ' open question. We suggest that the inflow is dominated by the relatively 

low angular momentum material originating in one of the nearby group 
of stars. Such material flows in directly towards the black hole up to the 
distance of order of ten Schwarzschild radii or less, where it hits the angular 
momentum barrier which leads naturally to a flow variability. We study 
both the analytical and the numerical solutions for the flow dynamics, and 
we analyze the radiation spectra in both cases using the Monte Carlo code 
to simulate the synchrotron, bremsstrahlung and the Compton scattering. 
Our model roughly reproduces the broad band spectrum of Sgr A* and 
its variability if we allow for a small fraction of energy to be converted 
^v^j ■ to non-thermal population of electrons. It is also consistent (for a range of 

viewing angles) with the strong constraints on the amount of circumnuclear 
material imposed by the measurements of the Faraday rotation. 
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1 INTRODUCTION 

The character of the accretion flow onto a black hole depends on the initial angu- 
lar momentum of the material. This angular momentum is specified by the outer 
boundary conditions which depend on the relative motion of the donor with respect 
to the black hole. This angular momentum corresponds to a certain circularization 
radius, i.e. the radius where this angular momentum is equal to the local Keplerian 
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NO-ISBN-SET-X © ???? - SU in Opava. All rights reserved. 



2 B. Czerny, M. Moscibrodzka, D. Proga, T. Das, and A. Siemiginowska 

value. In binary systems the material comes from the secondary star and in general 
is possess high angular momentum due to the orbital motion. In Low mass X-ray 
binaries the flow proceeds through an inner Lagrange point and the circularization 
radius is a significant fraction of a Roche radius around a black hole, of order of 
10 4 i? g (R g = GM/c). In high mass X-ray binaries the accretion flow comes from 
the intercepted focused wind, so the circularization radius is smaller but still large, 
of order of 10 3 i? g . In such case the inflowing material form an accretion disk around 
a black hole and the inflow proceeds due to the angular momentum transfer. Apart 
from the outermost region and the region close to the ISCO (innermost stable cir- 
cular orbit), the distribution of the angular momentum is relatively smooth and not 
much different from the Keplerian law. The exact departures from Kcplcrian motion 
depends on the disk temperature (or more exactly, on the pressure distribution). 

In active galactic nuclei (AGN) the source of material is less specified. The 
material comes either from the stars (in the form of stellar winds) or from the 
gaseous phase of the galactic material. Bright AGN (quasars, Seyfert 1 galaxies) 
show the presence of accretion disks similar to the disks in binary systems so we can 
conclude that the angular momentum reaching the galactic center is high. In sources 
showing water maser activity we observe the outer parts of the disk directly, and 
in most sources the motion of the disk material is Keplerian. However, in weakly 
active galaxies like Sgr A* or giant elliptical galaxies we see no direct evidence of 
a disk. In Sgr A* the presence of the cold disk is actually excluded by the lack 
of eclipses of the stars which move very close to the central black hole and are 
systematically monitored since several years. 

Since in weakly active galaxies there are no direct observational arguments for 
any value of the angular momentum of the donated material and the location of 
material sources, three types of models are being considered: 

• high angular momentum flow, with circularization radius of order of hundreds- 
thousands of Rg 

• low angular momentum flow, with circularization radius of order of a few R g 

• spherical and quasi-spherical accretion, without angular momentum barrier. 

The high angular momentum flow solutions for weakly active galaxies gener- 
ally belong to ADAF (advection dominated accretion flow) family (Ichimaru, 1977; 
Narayan and Yi, 1994), with possibly additional effects like outflows (Blandford and 
Begelman, 1999) and convection. In this case the flow is not exactly Keplerian since 
the pressure gradients are important, but the local ratio of the angular momentum 
to the Keplerian angular momentum in most part of the flow is not wildly different 
from unity, and the angular momentum transfer (through viscosity) or angular mo- 
mentum loss (through magnetic wind) at all radii is essential. Stationary solutions 
usually exist, and asymptotically the density of the flow approaches zero at infinity. 

In spherical and quasi-spherical flow there is no angular momentum barrier so 
the loss of angular momentum is not the necessary condition for the accretion to 
occur. Examples of such solutions are: purely spherical Bondi flow or flows where 
the angular momentum density is below the minimum angular momentum at the 
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circular orbit around a black hole which is given by 

Imin = 3V3GM/C (1) 

in case of Schwarzschild black hole; more general formula for a Kerr black hole 
can be found in (Bardeen et al., 1972). In Bondi solution (Bondi, 1952; Baganoff 
et al., 2003) the outer boundary condition arc specified by the density and the 
temperature of the uniform medium surrounding black hole at large distances. The 
flow velocity is zero at infinity, the inflow becomes transonic at the Bondi radius, 
and the supersonic flow reaches the black hole horizon. The Bondi radius depends 
significantly on the gas properties (e.g. politropic index; Bondi radius is of order of 
thousands of R g for relativistic flow with 7 = 4/3 but is approaches zero if 7 — > 5/3, 
typical for perfect fluid non- relativistic solution), but the accretion rate is much less 
sensitive to those assumptions. Purely Bondi flow has generally very low radiative 
efficiency so it cannot reproduce the observed luminosity in most weakly active 
galaxies (Moscibrodzka, 2006). If the accreting material at the outer boundary 
condition has certain angular momentum I < l m i n , the dynamics of the flow is 
slightly modified in comparison with Bondi flow and the flow is not spherically 
symmetric any more but the stationary solution for the flow always exists. 

The intermediate case of low angular momentum the situation is the most com- 
plex as initially the flow behaves as the Bondi flow but close to the black hole the 
flow starts suddenly to feel the angular momentum barrier (Abramowicz and Zurek, 
1981). In this case analytical stationary solutions frequently do not exist. In nu- 
merical solutions the flow is variable and does not reach a stationary solution in the 
computing time. If the angular momentum of the donated material is also a subject 
of changes (e.g. the result of the stellar motion) , a truly stationary solution indeed 
can never be reached for physical reasons. 

In the case of Sgr A* the available spatial resolution is the highest and we can have 
the best insight into the sources of material (Genzel and Karas, 2007). Therefore, 
in the present paper we concentrate specifically on this source and we argue that 
the low angular momentum flow is an interesting and promising option for the flow 
description. 



2 SGR A* SURROUNDING AND THE SOURCE OF MATERIAL 

At the central parts of our Galaxy there are a few stellar populations and each of 
them provides stellar winds. The closest O/B stars (at distances of a small fraction 
of a parsec) used for the mass measurements have moderate winds; much stronger 
winds come from more evolved stars being at distances up to a few parsecs. If 
a single star dominates as a donor star, and the wind velocity is larger than the 
stellar orbital velocity, there is a zero angular momentum line joining the donor 
star and the black hole, so the net angular momentum flowing in is likely to be low 
(Loeb, 2004; Moscibrodzka et al., 2006). Estimates of the ram pressure indicate 
that a single source, IRS 13 E (the compact group of Wolf-Rayet stars) with the 
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strongest wind (e.g. Rockefeller et al. (2004)) indeed dominates independently from 
the relatively large distance from the central black hole (Moscibrodzka et al., 2006). 

Although the ram pressure argument strongly depends on the adopted wind out- 
flow rate and wind velocity, the additional argument against the dominance of the 
nearby young O/B stars comes from the lack of obvious correlation between the 
activity level and stellar passages. On the other hand, un-modulatcd flow can come 
from stars forming a mysterious ring-like structure at a distance of a parsec scale 
(Paumard et al., 2006). In this case the material might have very large angular 
momentum, and even a cold disk may form, as postulated by Nayakshin and Sun- 
yaev (2003). However, the periodicity seemed to be seen in the NIR and X-ray flares 
(Genzel et al., 2003; Belanger et al., 2006) and the absence of eclipses does not seem 
to support cold disk scenario. The motion of IRS 13E differs from that of other 
stars and seems to be significantly eccentric. Detailed studies of the observational 
consequences of all scenarios are necessary to solve the issue. 



3 ANALYTICAL SOLUTIONS FOR THE DYNAMICS 

Analytical solutions can be obtained if the description of the flow is simplified: 
outflow and viscosity are neglected, the flow is politropic and the gravity is described 
by pseudo-Newtonian potential. The problem was studied in numerous papers 
following the idea of Abramowicz and Zurek (1981). 

3.1 Transonic solutions with and without shocks 

The existence and the character of solution depends critically on the interplay of 
the solution parameters: politropic index, energy density and angular momentum 
density, e.g. Das et al. (2003). The equations usually show the presence of three 
critical points: the outer one is direct generalization of the Bondi radius, the inner 
one describes the expected transition to the supersonic flow below the ISCO and 
the intermediate point is unphysical (in a sense it is not of the saddle type and the 
flow cannot pass there smoothly from subsonic to supersonic solution). Usually the 
flow passes from subsonic flow at infinity to supersonic flow close to horizon either 
at the outer or at the inner critical point. For a certain range of the parameters, 
it is also possible to find a second solution, with a shock (satisfying the standard 
Rankine-Hugoniot conditions) located between the outer and the inner sonic point. 
Such a shock may, but not must, form. 

The eventual shock development is likely to be related to the past state of the 
accretion flow. For a fixed value of politropic index and energy, there is a specific 
value of the angular momentum density at which the solution changes from transonic 
at outer to transonic at inner critical points. This change of flow properties is 
dramatic, and therefore the flow slowly crossing this angular momentum border is 
likely to develop a shock instead of following a new shock-less solution since this 
allows for a slow and continuous change in the flow properties (shock is initially 
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Figure 1. The two exemplary spectra models of SGR A* for two values of the angular 
momentum density in analytical model. One of the solutions extends only down to 20R g 
and correspond naturally to the weakly variable part of the flow. Thermal distribution of 
electrons was assumed. Dots mark the data points or upper limits in radio and NIR, and 
three representative levels of radio emission from Chandra are marked as power laws with 
slope errors. 



weak). The discussion of this issue will be presented elsewhere (Das, Czerny & 
Moscibrodzka, in preparation). 

For large enough angular momentum there is no analytical solution which extends 
down to the black hole even if a shock is allowed, and in this case the dynamics can 
be described analytically only e.g. down to a few a a few tens of R g (Moscibrodzka 
et al., 2006). This part of the flow is likely to be stationary while the inner part 
must form a kind of unstable ring. 



3.2 Exemplary spectra 

The radiation spectra were calculated with the code described in Moscibrodzka 
(2006) and later generalized to non-spherical distribution of the inflowing material 
(Moscibrodzka et al., 2006). The emission processes taken into account include 
synchrotron radiation, bremsstrahlung and Compton scattering. 

Two examples of the spectra are shown in Fig. 1. The spectra reproduce the 
NIR peak although are short of the data points at long wavelengths. The level of 
the continuous emission in X-rays is reproduced by the higher angular momentum 
model, so we can speculate that the emission of the innermost unstable ring, not 
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described by a stationary solution, may account for the strongly variable part of 
the X-ray flux. 



4 MHD SIMULATIONS OF THE FLOW DYNAMICS 

4.1 Flow description 

Exploratory phase of the use of MHD simulations for modeling Sgr A* time-dependent 
spectra made use of the simulations which were performed by Proga and Bcgclman 
(2003). The setup assumed the almost Bondi flow. The input of material was set 
at 1.2 Rb (where Rb is the Bondi radius), which in turn was equal to 2 x 10 3 R g , 
and the angular momentum of the new material was AR g c at the equatorial plane. 
The time-dependent computations were performed for a few dynamical timescales 
at the outer radius which gave the time-dependent distribution of the energy den- 
sity, velocity and the magnetic field. The flow was strongly time-dependent, with 
large scale fluctuations seen at the end of the computer run as well. The inflow was 
accompanied by a significant outflow. The computations never fully reached the 
exact stationarity condition in a sense that the inflow rate was not quite balanced 
by the outflow and accretion rate. However, in reality the flow also may not have 
time to reach equilibrium since the donor star moves (long term effects) and the 
Wolf-Rayet stellar winds show considerable dumpiness and overall variability (short 
timcscalc effects). 

Time-dependent electron temperature distribution was obtained by assuming sev- 
eral channels of energy transmission: (i) Coulomb coupling between the ions and 
electrons, (ii) direct heating of electrons due to compression (iii) electron radiative 
cooling (iv) electron advection. Optionally, we also allowed a fraction of energy to 
be in a form of non-thermal electrons. In order to simplify the computations, at 
each moment the flow was assumed to be stationary i.e. time-dependent solutions 
for the flow dynamics were treated as frozen frames. For a detailed description, see 
Moscibrodzka et al. (2007). 

4.2 Spectral variability 

Flow variability was reflected in variability of the broad band spectra. Exemplary 
states are shown in Fig. 2. We sec that the pure thermal electron distribution cannot 
represent the X-ray variability. The variable emission comes from the inner region 
of the flow in the form of synchrotron emission, and for thermal electrons this 
component does not extend to X-rays. Brcmsstrahlung emission, and Compton 
scattered emission comes from more extended region where variability is weaker 
and/or smeared. Therefore, the presence of non-thermal population is essential. 

The amplitude of the variability in the timescales recorded in MHD simulations (a 
day) was very large, over an order of magnitude so the variability in the overall ac- 
cretion rate (one order of magnitude) is additionally enhanced at some wavelengths 
(particularly at NIR) by the spectral effects. 
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Figure 2. MHD simulations: the exemplary broad band spectra of SGR A* for four 
different time moments labeled as A,B,C and D. Dots mark the data points or upper 
limits in radio and NIR, and three representative levels of radio emission from Chandra 
are marked as power laws with slope errors. Lower lines represent solutions with only 
thermal electrons, upper lines show solutions with a fraction of energy in a population of 
non-thermal electrons. 



The relation between the dynamics and the spectra is rather complex as it is 
strongly non-linear. There is no one-to-one correspondence between the accretion 
state and the predicted spectrum. For example, two states with similar accretion 
rates at the inner edge can have different spectra, or two states with different 
accretion rates can have similar spectra. 

The specific issue is the description of the outer parts of the flow. Bremsstrahlung 
dominates there, and the measured flux depends on the size of the emitting region. 
Since the spatial resolution available in X-ray band is low in comparison with the 
modeled region (even in Sgr A*!) this is an additional problem in comparing the 
models to the data. 



4.3 Faraday rotation measure 

The strongest observational constraint for the models comes from the estimates 
of the Faraday rotation measure towards Sgr A*. The observed polarization and 
the change of its angle indicate very low density plasma along the line of sight to 
the source. Our MHD model allows to calculate the integrated expected rotation 
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Figure 3. The determination of the Faraday rotation measure from the model, as a 
function of the viewing angle of an observer, for the same four different time moments 
labeled as A,B,C and D. Horizontal lines show the upper and the lower limit from the 
measurements of Marrone et al. (2006). 



measure as a function of the line of sight. The result is shown in Fig. 3. For some 
inclination angles the model is consistent with the observational limits. 



5 RECENT DEVELOPMENTS 
5.1 GR effects 

The computations presented in previous sections did not include the effects of Gen- 
eral Relativity properly — either purely Newtonian approach was used (with GR 
effects mimicked by the flow cut-off at 6R g ) or pseudo-Newtonian potential was 
adopted. Since the black hole in Sgr A* is likely to be rotating, more appropriate 
approach would be useful, or at least some estimates must be performed of an error 
to the spectrum due to the negligence of GR effects. 

Analytical solutions for the low angular momentum flow in the Kerr metric are 
well known (Barai et al., 2004) and can be used to determine the flow dynamics. 
Spectral computations using the Monte Carlo code require computing millions of 
photons so ray tracing for each of them is too time consuming if a number of models 
is to be calculated. However, we can calculate the escape probability of a photon as 
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Figure 4. The radial dependence of the Mach number for a solution with the shock and 
the corresponding escape probability; Kerr parameter a = 0.3. Dashed line in the lower 
panel shows the escape probability for a solution without a shock, with much larger radial 
velocity of the emitter in the innermost part of the flow. 



a function of radius, for a given Kerr parameter and we can test whether it depends 
significantly on the motion of the emitter (i.e. the solution for the flow dynamics). 

The computations are done in a standard way (see e.g. Schee et al. (2005) and 
the references therein) by integrating the photon paths for 100 000 photons emitted 
isotropically within the frame of the emitter at a given radius, and the calculations 
are performed for several radial points. 



10 B. Czerny, M. Moscibrodzka, D. Proga, T. Das, and A. Siemiginowska 



This work is still in progress, but in Fig. 4 we show an example for a specific 
dynamical solution with a shock. For a comparison, in the lower panel we also show 
the result for the solution without a shock (dashed line). We see that the material 
effective emissivity is much higher in the case with shock, even without the shock 
emission included. Therefore, the GR effects should be estimated separately for 
each dynamical solution. 

5.2 New dynamical simulations 

Also the dynamical MHD simulations used in the previous study was not fully 
satisfactory from the point of view of Sgr A* modeling. First of all, the results of 
the simulations were not recorded densely enough in time to allow us to follow as 
fast time variability as 17 minutes seen in QPO-likc events. New simulations, better 
suited for the center of our Galaxy, are currently being performed (Moscibrodzka 
& Proga, in preparation). The preliminary results show that indeed the shortest 
timescale variability is seen in the dynamical simulations. 



6 CONCLUSIONS 

Low angular momentum accretion flow is a promising scenario for the accretion onto 
Sgr A* due to its natural variability pattern. The flow is slightly more energetically 
efficient than the purely spherical Bondi flow and can reproduce both the required 
level of the luminosity and is consistent with the data on Faraday rotation measure. 
The overall broad band spectra are also roughly reproduced if a fraction of energy 
is allowed to be converted the non-thermal population of electrons. The current 
results are therefore encouraging, and the further work is in progress. 



ACKNOWLEDGEMENTS 

The present work was supported by the Polish Grant 1P03D 008 29 and the Polish 
Astroparticle Network 621/E-78/SN-0068/2007. 



REFERENCES 

Abramowicz, M. A. and Zurek, W. H. (1981), Rotation-induced bistability of transonic 

accretion onto a black hole, ApJ, 246, pp. 314-320. 
Baganoff, F. K., Maeda, Y., Morris, M., Bautz, M. W., Brandt, W. N., Cui, W., Doty, 

J. P., Feigelson, E. D., Garmire, G. P., Pravdo, S. H., Ricker, G. R. and Townsley, L. K. 

(2003), Chandra X-Ray Spectroscopic Imaging of Sagittarius A* and the Central Parsec 

of the Galaxy, ApJ, 591, pp. 891-915, arXiv:astro-ph/0102151. 
Barai, P., Das, T. K. and Wiita, P. J. (2004), The Dependence of General Relativistic 

Accretion on Black Hole Spin, ApJL, 613, pp. L49-L52, arXiv:astro-ph/0408170. 



Spinning test particles in equatorial plane of KdS spacetimes 11 



Bardeen, J. M., Press, W. H. and Teukolsky, S. A. (1972), Rotating Black Holes: Locally 
Nonrotating Frames, Energy Extraction, and Scalar Synchrotron Radiation, ApJ, 178, 
pp. 347-370. 

Belanger, G., Terrier, R., de Jager, O. C, Goldwurm, A. and Melia, F. (2006), Periodic 

Modulations in an X-ray Flare from Sagittarius A*, Journal of Physics Conference 

Series, 54, pp. 420-426, arXiv : astro-ph/0604337. 
Blandford, R. D. and Begelman, M. C. (1999), On the fate of gas accreting at a low rate 

on to a black hole, MNRAS, 303, pp. L1-L5, arXiv : astro-ph/9809083. 
Bondi, H. (1952), On spherically symmetrical accretion, MNRAS, 112, pp. 195 — h- 
Das, T. K., Pendharkar, J. K. and Mitra, S. (2003), Multitransonic Black Hole Accretion 

Disks with Isothermal Standing Shocks, ApJ, 592, pp. 1078-1088, arXiv : astro-ph/ 

0301189. 

Genzel, R. and Karas, V. (2007), The Galactic Center, in V. Karas and G. Matt, editors, 

IAU Symposium, volume 238 of IAU Symposium, pp. 173-180. 
Genzel, R., Schodel, R., Ott, T., Eckart, A., Alexander, T., Lacombe, F., Rouan, D. and 

Aschenbach, B. (2003), Near-infrared flares from accreting gas around the supermassive 

black hole at the Galactic Centre, Nature, 425, pp. 934-937, arXiv : astro-ph/0310821. 
Ichimaru, S. (1977), Bimodal behavior of accretion disks - Theory and application to 

Cygnus X-l transitions, ApJ, 214, pp. 840-855. 
Loeb, A. (2004), Direct feeding of the black hole at the Galactic Centre with radial 

gas streams from close-in stellar winds, MNRAS, 350, pp. 725-728, arXiv : astro-ph/ 

0311512. 

Marrone, D. P., Moran, J. M., Zhao, J.-H. and Rao, R. (2006), Interferometric Measure- 
ments of Variable 340 GHz Linear Polarization in Sagittarius A*, ApJ, 640, pp. 308-318, 
arXiv:astro-ph/0511653. 

Moscibrodzka, M. (2006), Spherical accretion in nearby weakly active galaxies, A&A, 450, 
pp. 93-103, arXiv:astro-ph/0512527. 

Moscibrodzka, M., Das, T. K. and Czerny, B. (2006), The pattern of accretion flow on to 
Sgr A*, MNRAS, 370, pp. 219-228, arXiv:astro-ph/0604516. 

Moscibrodzka, M., Proga, D., Czerny, B. and Siemiginowska, A. (2007), Accretion of Low 
Angular Momentum Material onto Black Holes: Radiation Properties of Axisymmetric 
MHD Flows, ArXtv e-prints, 707, 0707.1403. 

Narayan, R. and Yi, I. (1994), Advection-dominated accretion: A self-similar solution, 
ApJL, 428, pp. L13-L16, arXiv:astro-ph/9403052. 

Nayakshin, S. and Sunyaev, R. (2003), Close stars and an inactive accretion disc in Sgr 
A*: eclipses and flares, MNRAS, 343, pp. L15-L19, arXiv:astro-ph/0302084. 

Paumard, T., Genzel, R., Martins, F., Nayakshin, S., Beloborodov, A. M., Levin, Y., 
Trippe, S., Eisenhauer, F., Ott, T., Gillessen, S., Abuter, R., Cuadra, J., Alexander, 
T. and Sternberg, A. (2006), The Two Young Star Disks in the Central Parsec of 
the Galaxy: Properties, Dynamics, and Formation, ApJ, 643, pp. 1011-1035, arXiv: 
astro-ph/0601268. 

Proga, D. and Begelman, M. C. (2003), Accretion of Low Angular Momentum Material 
onto Black Holes: Two-dimensional Magnetohydrodynamic Case, ApJ, 592, pp. 767- 
781, arXiv:astro-ph/0303093. 

Rockefeller, G., Fryer, C. L., Melia, F. and Warren, M. S. (2004), Diffuse X-Rays from the 
Inner 3 Parsecs of the Galaxy, ApJ, 604, pp. 662-670, arXiv:astro-ph/0309497. 

Schee, J., Stuchlfk, Z. and Juran, J. (2005), Light escape cones and raytracing in Kerr 



12 B. Czerny, M. Moscibrodzka, D. Proga, T. Das, and A. Siemiginowska 



geometry, in S. Hledfk and Z. Stuchlfk, editors, RAGtime 6/7: Workshops on black 
holes and neutron stars, pp. 143-155. 



